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ABSTRACT 

We combine time-dependent multi-waveband flux and linear polarization observations with sub- 
milliarcsecond-scale polarimetric images at A = 7 mm of the BL Lacertae-type blazar OJ287 to locate the 
7-ray emission in prominent flares in the jet of the source > 14 pc from the central engine. We demonstrate 
a highly significant correlation between the strongest 7-ray and millimeter- wave flares through Monte-Carlo 
simulations. The two reported 7-ray peaks occurred near the beginning of two major mm- wave outbursts, each 
of which is associated with a linear polarization maximum at millimeter wavelengths. Our Very Long Baseline 
Array observations indicate that the two mm- wave flares originated in the second of two features in the jet that 
are separated by > 14 pc. The simultaneity of the peak of the higher- amplitude 7-ray flare and the maximum in 
polarization of the second jet feature implies that the 7-ray and mm- wave flares are co- spatial and occur > 14 pc 
from the central engine. We also associate two optical flares, accompanied by sharp polarization peaks, with the 
two 7-ray events. The multi- waveband behavior is most easily explained if the 7-ray s arise from synchrotron 
self-Compton scattering of optical photons from the flares. We propose that flares are triggered by interaction 
of moving plasma blobs with a standing shock. The 7-ray and optical emission is quenched by inverse Comp- 
ton losses as synchrotron photons from the newly shocked plasma cross the emission region. The mm-wave 
polarization is high at the onset of a flare, but decreases as the electrons emitting at these wavelengths penetrate 
less polarized regions. 

Subject headings: Galaxies: active — galaxies: jets — BL Lacertae objects: individual (OJ287) — radio 
continuum: galaxies — gamma rays: general — polarization 



1. INTRODUCTION 

The information that 7-ray observations can provide on 
the physical properties of active galactic nuclei in gen- 
eral, and blazars in particular, depends on where such 7- 
ray emission originates, which is still under debate (e.g., 
iMarscher & Jorstadll2010|) . Variability studies of 7-ray and 
lower frequency emission from blazars provide important in- 
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sights into this problem by relating time scales of variability 
with physical sizes in the source. However, although 7-ray 
variability in blaz ars can occur on very short time- scales (o f 
a few hours, e.g.. iMattox et alJ 119971: iFoschini et aPl2Q10h . 
this does not necessarily imply that high-energy flares take 
place at short distances (<g j 1 pc) fro m the central engine (see 
IMarscher & Jorstad 20 1 0). In f act, iLahteenmaki & Valtaojal 
(120031) and I Jorstad et al.1 (l2001a llbh provided evidence that 7- 
ray outbursts may arise in the mm-wave emitting regions > 
1 pc downstream of the central engine. The Large Area Tele- 
scope (LAT) onboard the Fermi Gamma-ray Space Telescope 
has sufficient sensitivity to test this by providing well- sampled 
light curves of blazars that allow detailed studies of the timing 
of 7-ray flares relative to those at other spectral ranges (e.g . , 
lAbdo et al.ll20l0dl: I Jorstad et al.ll20lol: IMarscher et alJl20TQh . 

The techn i que d eveloped by IMarscher et al. I (120101) and 
iJorstad etail (120101) uses ultra-high angular-resolution ( ~ 
0.15milliarsecond) monitoring with very long baseline inter- 
ferometry (VLBI) to resolve the innermost jet regions and 
monitor changes in jet structure. Observations with roughly 
monthly observations, supplemented by more concentrated 
campaigns, provide time sequences of total and polarized in- 
tensity images of the parsec- scale jet that can be related to 
variations of the flux and polarization at higher frequencies. 

In this paper, we employ this technique to investigate 
the location of the flaring 7-ray emission in the BL Lacer- 
tae (BL Lac) object OJ287 (z = 0.306), a well studied and 
highly variable bl azar at all available spectral ranges (e.g., 
lAbdo et aTl l2010alcT) . Throughout this paper we adopt the 
standard ACDM cosmology with Hq=71 km s -1 Mpc -1 , Qm = 
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0.27, and Q\ = 0.73, so that 1 milliarcsecond (mas) corre- 
sponds to a projected distance of 4.48 pc, and a proper motion 
of 1 mas/yr corresponds to a superluminal speed of 19 c. 

2. OBSERVATIONS 

Our polarimetric observations of OJ287 include (1) 7 mm 
(43 GHz) Very Long Baseline Array (VLB A) images (Fig.Q]), 
mostly from the Boston University blazar monitoring pro- 
granO (2) 3 mm (86 GHz) monitoring with the IRAM 30 m 
Telescope, and (3) optical (R and V band) photo-polarimetric 
observations from several observatories (Figs. [2] and [3]). The 
optical facilities include Calar Alto (2.2 m telescope, obser- 
vations under the MAPCAlQ program), Steward (2.3 and 
1.54 m telescope JI 7 !), Lowell (1.83 m Perkins Telescope), St. 
Petersburg State University (0.4 m telescope), and Crimean 
Astrophysical (0.7 m telescope) observatories. To these we 
add R-band po larimetric (and V-band photometric) data from 
IVillforth et all d2010h . The total flux light curves ana- 
lyzed here (see Fig. [2]) are from the Fermi-LAI 7-ray (0.1- 
200 GeV) and Swift X-ray (0.3-10 keV) and optical (V-band) 
data available from the archives of these missions, by the Yale 
University SMARTS prograrrQ by the Submillimeter Array 
(SMA) at 1.3 mm (230 GHz) and 850 /im (350 GHz), by the 
IRAM 30 m Telescope at 1 .3 mm, and by the Metsahovi Radio 
Observatory 14 m Telescope at 8 mm (37 GHz). 

Our data analysis follo ws the procedures f rom previ- 
ous studies: (1) VLB A: IJorstad et al l (|2005l) ; (2) op- 
tical polarimetric d ata: IJorstad et al.l (120101); (3) IRAM 
data: lAgudo"et~aTl $2(MA20]^ ^ et all 

d2007h: (5) Metsah ovi: iTerasranta et al.l figp); (6) Swift 
IJorstad et al.l (120101) : and (7) Fermi LAT: iMarscher et al.1 
(120101) . For the LAT data, we considered a 15° radius cen- 
tered on OJ287 and used the maximum-likelihood routine GT- 
LIKE to model the 0.1-200 GeV spectra from this object, 
CRATES J0856+2057, OJ 290, and 1FGL J0902.4+2050 as 
single power laws. We measured the flux in 7 -day bins and 
fixed the slope of the photon spectrum of OJ287 at -2.4, (as 
determined over th e first 1 1 months of LAT observations by 
lAbdo et alJl2010bh . 

3. ANALYSIS OF OBSERVATIONAL RESULTS 

3.1. 7 mm Jet Structure and Kinematics 

Fig. Q] reveals a radically different sub-milliarcsecond- 
scale jet structure a t 7 mm than found previously in OJ287. 
IJorstad etaLl (120051) reported a bright core with a one-sided 
jet at a position angle PA w -110° in 1998-2001. Our 
new images show the innermost jet region oriented with PA 
from - -10° to - -40° between 2008 and 2010, among the 
most ext reme shifts in app arent jet direction yet observed 
(see e.g., Ag udo et aLll2007l) . The change in orientation con- 
flicts with the 12-yr period precession model presented by 
lTateyama~ & Kingham (2004) that predicted the jet to be at 
PA w -95° in 2009. A more detailed study of this issue and 
of the origin of jet wobbling in OJ287 will be presented else- 
w here. 

IJorstad etaLl (120051) also reported a set of quasi- stationary 
features at average locations ~ 0.1, ~ 0.3, and ~ 1 mas from 
the innermost jet feature in OJ287. Fast superluminal fea- 
tures propagating downstream with speeds > 10 c were also 

15 http : //web .bu . edu/blazars/VLBApro ject . html 

16 http : //www. iaa . es/~iagucio/research/MAPCAT 

17 Data listing: http : // james . as . arizona . edu/~psmith/Fermi 

18 Data listing: [Ft tp : / /www . astro . yale . eciu/sma rts/glast| 
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Fig. 1.— Sequence of 7 mm VLBA images of OJ287 in 2008-2010. Images 
are convolved with a FWHM = 0.15 mas circular Gaussian beam. Contour 
levels represent 0.2, 0.4, 0.8, 1.6, 3.2, 6.4, 12.8, 25.6, 51.2, 90.0% of the peak 
total intensity of 6.32Jy/beam. The color scale indicates linearly polarized 
intensity, whereas superimposed line segments represent the orientation of 
the polarization electric-vector position angle. 
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observed to move across these quasi- stationary features. We 
model the brightness distribution of the source at 7 mm with 
a small number of circular Gaussian components. Our model 
fits include a bright quasi- stationary feature (CI) ~ 0.2 mas 
from the innermost jet region (CO). The identification of CO 
as the innermost jet feature is justified by the decreasing inten- 
sity westward of CI, and by the detection after 2010 March of 
superluminal motion of features Ml and M2 toward the west- 
southwest of CI with speeds of 10.8c ± 1.3c and 6.7c ± 1.4c) 
and possibly M3, which crossed CI in 2010 Oct. (preliminary 
speed of > 10c). 

3.2. Flares in the CI Jet Region at 1 mm and 7 mm 

Before Ml, M2, and M3 appeared in the jet, CO and CI 
typically contained > 90% of the 7 mm flux, and therefore 
governed the mm- wave evolution, of OJ287. In particular, the 
two most prominent 1 mm flares ever reported in OJ287 (Amm 
and B mm as labeled in Fig. [2]) took place in CI, as indicated by 
the correspondence of events in the 7 mm light curve of this 
component with those at other millimeter wavelengths. The 
multiple peaks of flare B mm may be related to the passage of 
Ml, M2, and M3 through CI, although we cannot verify this 
given the uncertainties in the trajectories of these knots. No 
moving features related to A mm are apparent in the images, 
although the polarization westward of CI in the 2009-02-22 
and 2009-06-21 images suggests that weak knots were present 
at those epochs. 

From the a ngle of the jet axis t o the line of sight in O J287 
(1?9-4?1; see lJorstad et aLll2005l: IPushkarev et al.ll2009h . and 
the mean projected separation of CI from CO at the time of 
start of A mm and B mm (0.23 ± 0.01 mas), we estimate that CI 
is located > 14 pc downstream of CO, the innermost jet region 
detected in our images. The actual distance between the cen- 
tral engine in OJ287 and CI must be even greater if CO lies 
downstre am of the acceleration and collimation zone of the 
jet (ACZi lJorstad et al.ll2007l:lMarscher et al.ll2008L l2010h . 

3.3. 7 -Ray Flares 

The two most pronounced 7-ray flares take place during the 
initial rising phases of and B mm . This conforms with th e 
pattern found in the 1990s by lLahteenmaki & Valtaojal (120031) . 
The two 7-ray flares, peaking on 2008 October 4 (A 7 ) and 
2009 October 24 (B 7 ), reached 0. 1-200 GeV photon fluxes a 
factor of ~ 2 and ~ 5, respectively, higher than the quiescent 
7-ray level (~ 1 x 10~ 7 phot/cm 2 /s). 

The discrete correlation function (DCF; lEdelson & Krolikl 
Il988|) between the 7-ray and 1 mm long-term light curve 
(Fig. [4]) possesses a prominent peak at a time lag ~ -80 days 
(7-ray leading). We compute the significance of this result 
relative to stochastic variabi l ity thr ough Monte Carlo simula- 
tions a s in IChatterjee et al.l (l2008|) and iMa x-Moerbec ket al.l 
(120101) . Our method follows iTimmer & Koenigl (119951) for 
the simulation of N = 5000 different pairs of 7-ray and 1 mm 
light curves characterized by the mean and standard devia- 
tion of the observed light curves, and by power-law shaped 
power spectral densities (PSD oc l/f a ) with slopes in a grid 
a 1 = {1.0,1.5,2.0} and a lmm = {1.0, 1.5,2.0,2.5 3.0|. The 
range o f a 1 values includes the means found by lAbdo et ail 
(l2010cT) . a 7 = 1.5 and 1.7 for quasars and BL Lacs, respec- 
tiv ely. For a imm, we choos e va lues covering the range derive d 
bv lChatterjee et al.l (I2008I) and lHufnagel & Bregmanl (I1992I) . 

The results of our simulations show that the DCF peak at 
a lag of ~ -80 days is significant at 99.7 % confidence in all 
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Fig. 2. — Light curves of OJ287 from mm- wave to 7-ray frequencies. The 
vertical lines denote the times of peak 7-ray flux of A 7 and B 7 . The solid 
lines in the next- to-last pane l represent fits to the two major 1 mm flares by 
using the method of lChatterieeeTall l2Q08h . RJD = Julian Date -2400000.0. 

of our simulations. This confirms the correlation between B 7 
and B mm , the most luminous j-ray and 1 mm flares in our 
data. The correlation of the 7-ray and 8 mm light curves in 
Figgis of similarly high significance. 

The optical light curves (especially V-band; Fig. O show 
two sharp flux increases at essentially zero time lag from A 7 
and B 7 . In contrast, the sparser time coverage of the 0.3- 
10 keV light curve does not allow us to make an unambiguous 
connection between the 7-ray and X-ray flares. 

3.4. Variability of Linear Polarization 

CO and CI dominate the evolution of the linear polariza- 
tion p and electric vector position angle \ at 7 mm in OJ287 
(Fig. [3]). However, whereas pco never exceeds 10%, CI ex- 
hibits the two largest peaks in p ever observed in OJ287 at 
7 mm, p C \ « 14% on 2008 November 4, and p C \ « 22% on 
2009 October 16. The first maximum in pa follows the peak 
of A 7 by one month, while the second more pronounced po- 
larization event is already in progress when the 7-ray flux of 
flare B 7 rises to a level of ~ 4 x 10~ 7 phot cm -2 s~ l . This 
coincidence of the strongest j-ray outburst and exceptionally 
strong polarization in CI identifies this feature > 14 pc from 
the central engine as the site of the variable j-ray emission. 
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Fig. 3. — Optical and mm-wave linear polarization of OJ287 as a function 
of time. The optical data includes R-band and 5000-7000 A observations. 
The similarity of the integrated linear polarization at 7 mm and 3 mm also 
allows us to combine them. Vertical lines are as in Fig. [2] 



The optical polarization peaks at essentially the same time 
as pci at 7 mm during flare B 7 . During both A 7 and B 7 , 
p op t ~ 35%, which requires a well-ordered magnetic field. 
However, comparable optical polarization levels also occur at 
other times. The shorter time scale and larger amplitude of 
variability of optical polarization, as compared with those at 
millimeter wavelengths, is consi stent with frequency depen- 
dence in the turbulence model of M arscher & J orstad ( 2 010|) . 

The optical and mm-wave linear polarization position an- 
gle is stable at x ~ 160°-170° (-20° to -10°)— similar to 
the structural position angle of the inner jet — both near A 7 
and B 7 and throughout most of the monitoring period. The 
corresponding direction of the magnetic field is transverse to 
the direction between features CO and CI. The stability of 
field direction is only greatly altered in the optical by spo- 
radic short-term rotations of y by up to 1 80° whe n the polar- 
ization is relatively low (see I Villforth et al.ll2010l) . consistent 
with th e behavior expected if the magnetic fi eld becomes tur- 
bulent dJonesHl98alD'Arcangelo et al.ll2007l) . 

3.5. Low Probability of Chance Coincidences 

We note that the correlation peak in Fig. [4] is dominated by 
events B 7 and B mm . Despite the prominence of flare A mm , 
the weakness of its 7-ray counterpart (A 7 ) prevents the pair 
of events from being apparent in the DCF at their time lag of 
~ -120 days. However, there is statistical support for the hy- 
pothesis that flares A 7 and A mm are also physically related. 
The probability that a 7-ray outburst (approximated as instan- 
taneous relative to the mm-wave flares) occurs by chance dur- 
ing a ~ 120-day rise time of a 1 mm flare peaking at > 5.5 Jy 
is only 17%. This probability is inferred from the duration 
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Fig. 4. — Discrete cross-correlation function between the 7-ray and 1 mm 
light curves of OJ287 (red points). The dotted curves at positive (nega- 
tive) DCF values denote 99.7% confidence limits for correlation rather than 
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of the long-term 1 mm SMA light curv43 and the number of 
such flares. Hence, the probability that two 7-ray flares at 
random times occur by chance during the rising phase of two 
mm-wave flares is 3%, i.e., events A mm and B mm are associ- 
ated with A 7 and B 7 , respectively, at 97.0% confidence level. 
This rises to 99.2 % if, instead of A mm and B mm , the two high 
pa mm-wave peaks are considered. 

4. DISCUSSION AND CONCLUSIONS 

The two 0. 1-200 GeV flares in OJ287 allow us to assess 
the correspondence between 7-ray and lower- frequency vari- 
ations. We find that two kinds of events at millimeter wave- 
lengths are related to these 7-ray outbursts at high signifi- 
cance: (1) the early, rising phases of the two most luminous 
1 mm flares ever detected in this blazar, (A mm and B mm ); and 
(2) two sharp increases to unprecedented levels of linear po- 
larization (~ 14% and ~ 22%) in bright jet feature CI > 14 pc 
from the central engine. These events also coincide with: (3) 
two sharp optical flares; (4) two peaks in optical polarization 
of ~ 35%; and (5) the similarity of optical and mm-wave po- 
larization position angle both during and between the flares at 
X~ 160°-170°. 

The exceptionally high polarization of CI during 7-ray flare 
B 7 provides extremely strong evidence that the event oc- 
curred in CI. This has two important implications. First, 
given the distance of CI from the central engine, the 0.1- 
200 GeV flares must be produced by inverse Compton (IC) 
scattering rather than nuclear collisions. Second, the 7-ray 
IC emission arises from either the synchrotron self-Compton 
(SSC) process or IC scattering of infrare d radiation from a 
hot, dusty torus of size ~ 10 pc (IC/dust; iBlazeiowsk i et al.l 
120041: ISokolov & Marscheil 120051) . An SSC model is possi- 
ble given the low ratio of 7-ray to synchrotron luminosity be- 
tween 10 14 and 10 15 Hz (« 2) in OJ287, based on the fluxes 
we measure and a spectral index of -1.5 between 10 14 and 

19 |http : // smal . sma . hawaii . edu/ callist /callist . html| 
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Fig. 5. — Sketch of the proposed model for the multi- wavelength flaring 
behavior of OJ 287. 



10 15 Hz (IVillforth et al.l 1201 Oh . Inverse-Compton scattering 
of dust emission predicts that the optical and 7-ray emission 
should vary toget her, since the electron energies involved are 
similar (see, e.g. jMarscher et al1l2010|) . This is the case dur- 
ing the 7-ray flares. The general data are therefore consistent 
with both the SSC and the IC/dust models. 

The observed behavior agrees with a scenario in which the 
optical and 7-ray flares are produced in CI by particle ac- 
celeration in a moving blob when it crosses a standing shock 
(Fig.© well beyond the ACZ. lD'Arcangelo et al.l (l2QQ7h have 
associated the innermost mm- wave jet emission feature with 
a recollimation shock at the end of the ACZ. We identify CO 
as such a feature, and CI as a second r e-collimation shock , 
as seen in hydrody namical simulations (iGomez et al.ll 19971: 
lAgudo et a l. 2001) and often in V LBA images djorstad et al.l 
I2001all2005|) . iJorstad et al.l (1201 0|) proposed a similar system 
of three conical shocks to explain the fine structure in the in- 
nermost 7 mm jet of 3C 454.3 in a feature that would be iden- 
tified as the "core" at centimeter wavelengths. In OJ287, the 
higher flux of CI relative to CO could be caused by curvature 
in the jet increasing the Doppler beaming with distance from 
CO so that it is maximized at CI. 

A blob is either a shock or another disturbance that prop- 
agates down the jet with significantly higher relativistic elec- 
tron density and magnetic field than in the ambient flow. The 
typical weakness of the observed polarization outside the high 
polarization peak s favors a turbulent plasma over a shock. 
ICawthornel (120061) shows that, in this case, the side of the con- 
ical shock nearest to the line of sight can be highly polarized 
with x parallel to the jet axis. The remainder of the conical 



shock, farther from the line of sight, has much lower polariza- 
tion. Because of light-travel delays, we first see the blob pen- 
etrate the near side, and therefore observe a major increase in 
polarization. As the outburst develops, more of the emission 
comes from the low-polarization far side, which decreases p 
while the mm- wave flux density continues to increase. This is 
the pattern observed during both flares. 

The mm- wave and optical flares start at essentially the same 
time as the magnetic field and electron energies near the lead- 
ing edge of the blob become amplified as they pass the stand- 
ing shock front. The mm- wave flux outburst continues as the 
relatively low-energy electrons fill the shocked region. The 
optical and 7-ray emission, produced by higher-energy elec- 
trons that cannot travel far before suffering radiative energy 
losses, is co nfined closer to the shock front where particles are 
accelerated dMarscher & Geaj [T985). If synchrotron losses 
were the sole factor, the optical and 7-ray flux would reach 
a plateau once the shocked blob plasma fills the entire layer 
behind the shock, declining only when the upstream side of 
the blob passes the shock front. However, if SSC emission 
dominates, the radiative losses will increase as synchrotron 
photons from the flare reach electrons that scatter them to 
high energies. The radiative energy losses then increase un- 
til the loss rate reaches a factor of [1 + (Lssc/^synch)] (~ 3 in 
the case of OJ287) higher than the initial case of only syn- 
chrotron losses. This decreases the volume of emission as a 
function of time in the optical and 7-ray ranges. The time- 
scale for this decrease is the light-travel time across the shock 
times (1 +z)S~ l ~ 20(£/19) -1 OV0.054mas) days, where a is 
the cross-sectional angul ar size of the jet and the Doppler fac- 
tor (5) is that derived by lJorstad et al.l (|2005|) . This is similar 
to the time-scale of the 7-ray flux decline of flare B 7 (~ 3-4 
weeks, see Fig. [2]). 

Quenching of the flares by increasing SSC energy losses 
operates only if the 7-ray s are produced by the SSC process. 
Hence, the multi-frequency behavior of the flares is difficult 
to reproduce in the IC/dust model. Also, infrared emission 
from the dusty torus has not been detected thus far in BL Lacs 
such as OJ287, as far as the authors know. We thus favor the 
SSC mechanism. 
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